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Abstract 

We reexamine Chew's method for extracting the ttNN coupling constant 
from np differential cross section measurements. Values for this coupling 
are extracted below 350 MeV, in the potential model region, and up to 1 
GeV. The analyses to 1 GeV have utilized 55 data sets. We compare these 
results to those obtained via x 2 mapping techniques. We find that these 
two methods give consistent results which are in agreement with previous 
Nijmegen determinations. 
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I. INTRODUCTION 



While the ttNN coupling constant has been known within a few percent for decades, 
there have been numerous attempts to more precisely determine its value, and to search 
for possible charge-independence breaking. Recent results from the Nijmegen NN and 
NN analyses jl], and the Virginia Tech ttN and iViV analyses 0, have found values 
of (7 2 / 47r near 13.7 with little evidence for charge dependence. However, the importance of 
this coupling continues to motivate independent studies |4| utilizing different methods of 
extraction. 

In 1958, Chew suggested || that the influence of the pion-pole might be measureable 
in NN elastic scattering, if appropriate kinematic regions were examined. This technique 
has recently been applied to pp — > nn and np elastic [0 differential cross section mea- 
surements. In these studies, the extracted values for g 2 /4:TT (g 2 /Air = 12.83 ± 0.36 || and 
g 2 /4:7T = 12.47 ± 0.36 0]) were somewhat below those found in Refs. [0], 0. 

In the present work, selected np differential cross section data Jt| — [15] below 1 GeV 



have been analyzed using a modified form of the Chew method f|, [11|. In the region below 



350 MeV, the "difference method" [|Tl|] was used in conjunction with the Nijmegen, Bonn 
and Paris potentials, and various partial-wave solutions. We considered differences between 
"model" and "experimental" values of the dimensionless quantity y = x 2 s (da / dQ) / (he) 2 
expanded in terms of a function of x = 1 + MT(1 + cos9)/p 2 , 

N 

V&t = Z/model + ^2ajbj(x) (1) 
j=0 

where s is the usual Mandelstam variable, T is the laboratory kinetic energy, and fi (M) is 
the pion (nucleon) mass. While it is conventional to expand (ym — y m odei) as a power series 
in x, we chose to expand in terms of Legendre polynomials, Pj(p), 

b = 1, bj(j >0) = Pj(p)-l, (2) 

with p = 1 — 2(x/x max ), which resulted in improved numerics but had no effect on the final 
results (several different expansions were made for comparison purposes). The difference 



between the model coupling constant and the experimental value is then given in terms of 
the first coefficient 

fexp = flmodel + °0> (3) 

in the above expansion. 

The results of our analyses are given in Tables I and II. In Table I, we list values of 
g 2 /47r found via the difference method applied to data below 1 GeV, using our partial- wave 
solution SM95 as input. In order to treat the different datasets as consistently as possible, 
we put a number of constraints on the fits. 

• Extractions were performed on data backward of 110 degrees. 

• Only those datasets with more than 15 qualifying points were accepted. 

• An upper limit for the variable x was taken to be 10. 

• The number of fitted parameters was given a lower limit of 4. 

The search was terminated either at a minimum x 2 / (degree of freedom) or when the error 
on g 2 /4ir exceeded 1. 

In order to gauge the effect of renormalizing datasets, two independent determinations 
of g 2 /4:7i were made for each experiment. In the first case, data were fitted without any 
renormalization. In the second case, data were first renormalized according to the input 
solution SM95. The average values of g 2 /4:7i, and their associated uncertainties, are listed in 
Table II for a number of input solutions and energy ranges. The uncertainties were calculated 
without accounting for systematic errors. We estimate that more realistic uncertainties are 
larger by at least a factor of 2. 

The values in Table I show considerable scatter, as is evident from Fig. 1. Given this vari- 
ability, the calculated averages are remarkably consistent. Coupling constants determined 
using potential model and partial-wave inputs below 400 MeV are in good agreement those 
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found using our partial-wave solutions up to 1 GeV. These averages are also quite consistent 
with the results of Refs. 0, §. 

As a check, we have repeated the coupling constant extractions of Refs. [[7|], (Tl[ and 
agree with their results. In Ref. [|11|] it is claimed that the 96 MeV and 162 MeV Uppsala 
cross sections ||, |TT[ imply a consistent value for g 2 /Air. This is not evident from the values 
quoted in Table I. We can, however, explain the discrepancy. While order 4 and order 6 fits 
to the data of Ref. || gave values of g 2 /Air consistent with those found in Ref. [|ll|], the order 



5 fit gave a somewhat larger value. In choosing the results listed in Table I, the order 6 fit 
was rejected, as the associated error was too large, according to our criterea for terminating 



a search. In contrast, the coupling derived from the data of Ref. |TT| was much less sensitive 
to the chosen order. 

As a further consistency check we mapped g 2 /An against \ 2 over different energy regions, 
analyzing pp and np data both separately and in combined fits. This extended a study 
reported earlier 0. In analyzing the combined dataset (pp and np) below 400 MeV, we found 
a value of coupling near 13.8. Results were obtained using the Nijmegen Coulomb distorted- 
wave Born approximation |^TJ for the one-pion exchange (OPE) contribution. A significantly 
higher value of g 2 /Air (near 14.4) was found if we used a plane-wave Born approximation 



for the OPE contribution. While this appears to confirm our Chew-extrapolation results, less 
consistency was found when the extraction was broken into separate pp and np contributions, 
or when the energy range was increased to 1.6 GeV. In fact the sign of (g 2 p — g 2 p ) changed 
when the analysis was extended from 400 MeV to 1.6 GeV. These results are listed in 
Table III. Results of our analyses to 1.6 GeV are displayed in Fig. 2. The \ 2 map from our 
analyses of elastic pion-nucleon scattering data is included for comparison. 

In summary, we have found that the value of g 2 /An extracted from NN scattering data is 
quite consistent with our previous determinations. The Nijmegen group has also generated 
X 2 maps [|TJ of the Bonner et al. |1(]] data below 350 MeV, obtaining results in qualitative 
agreement with ours. However, given the variation of g 2 /An values displayed in Fig. 1, we 
cannot claim to have improved upon the results found from our application of dispersion 



relations to elastic pion-nucleon scattering data. 
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Figure captions 



Figure 1. Plot of (? 2 /47r determinations from np differential cross sections versus the 
laboratory kinetic energy. The input solution was SM95. The solid line gives the 
average value. Results for unnormalized data are displayed. 

Figure 2. x 2 maps from fits to all NN data below 1.6 GeV (dot-dashed line), pp data 
only (dashed line), and np data only (dotted line). The result from fits to elastic 
pion-nucleon scattering data (solid line) is display for comparison. The variation 
of x 2 from minimal values (Sx 2 ) is plotted against g 2 j^n. 
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Table I. Values of g 2 /^n extracted using the difference method applied to data below 
1 GeV at backward angles (see text). The laboratory kinetic energy (T), the angular range 
(9) and number (N) of data used in the analysis, the order (M) of the fit, and the corre- 
sponding x 2 /(degree of freedom) are given for each dataset. In the first column of g 2 /Ait 
values, datasets were fit without renormalization. In the next column, a normalization factor 
(Norm), determined from solution SM95, was applied prior to fitting. 



T 

l^lvlt: V J 




a 

u 

\ ) 


N Ref 

M 


X l ao J 


(unnormalized) 

g 2 /47r Norm 


(normalized" 


91.0 


119 


- 177 


16 






4 


0.81 


12.59 


'1.41) 


1.056 


10.96 


'2.36) 


96.0 


117 


- 179 


32 


d 




4 


0.61 


16.50 


0.57) 


0.982 


16.37 


0.56 


162.0 


122 


- 178 


43 


1 


u 


4 


1.51 


12.96 


0.75 


1.050 


13.19 


0.77) 


162.0 


119 


- 179 


31 


11 


4 


1.05 


14.65 


;o.36) 


0.938 


14.32 


0.34 


177.9 


122 


- 179 


44 


10 


4 


1.08 


13.06 


;o.67) 


1.042 


13.25 


'0.69 


194.5 


121 


- 179 


42 


10 


4 


1.54 


12.15 


;o.63) 


1.040 


12.29 


0.65 


211.5 


120 


- 178 


43 


10 


4 


0.76 


12.75 


;o.54; 


1.029 


12.86 


0.55 


212.0 


111 


- 177 


30 


12 


6 


0.69 


14.67 


,0.68, 


0.970 


13.88 


0.49 


229.1 


120 


- 178 


49 


10 


4 


1.36 


12.80 


0.47 


1.027 


12.90 


0.48) 


247.2 


119 


- 178 


53 


10 


4 


0.73 


13.21 


0.42 


1.014 


13.26 


> 0.42) 


265.8 


118 


- 179 


63 


m 


4 


0.84 


12.97 


0.37 


1.002 


12.98 


0.37 


284.8 


117 


- 179 


73 


10 


4 


1.14 


13.32 


0.23 


1.024 


13.41 


0.24) 


304.2 


116 


- 179 


79 


10 


6 


0.97 


13.16 


0.76 


0.975 


12.72 


} 0.88 


319.2 


117 


- 177 


40 


12 


6 


0.85 


11.85 


0.69 


0.989 


11.75 


0.69 


324.1 


115 


- 179 


80 


10 


4 


1.10 


13.57 


0.18 


1.028 


13.67 


0.19 


344.3 


118 


- 179 


74 


10 


4 


1.12 


13.14 


0.21 


1.003 


13.15 


0.21 


365.0 


120 


- 179 


62 


10 


5 


1.33 


12.72 


;o.4e; 


1.046 


12.89 


0.50 


386.0 


121 


- 179 


55 


10 


4 


1.44 


13.51 


;o.2o) 


1.066 


13.73 


;o.2i) 
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Table I (continued). 
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Table II. Table of average g 2 /^ values obtained using partial-wave analyses (SM95, 
VV40, and NI93), and potential models (NY93, BONN, and PARIS) as input. N is the 
number of datasets involved in the analysis. 



Solution 


(Model) 
5 2 /4tt 


Ref 


Range 
(MeV) 


N 


(unnormalized) 


(normalized) 


SM95 


13.75 


present 


- 1000 


55 


13.58 (0.04) 


13.75 (0.04) 


SM95 


13.75 


present 


0-400 


18 


13.47 0.08 


13.51 0.08 


VV40 


13.75 


present 


0-400 


18 


13.55 0.08 


13.66 0.08 


NI93 


13.58 


[17| 


0-350 


16 


13.51 (0.08) 


13.65 (0.09) 



NY93 


13.58 


18 





-350 


16 


13.04 


BONN 


14.40 


19 





-325 


15 


13.72 


PARIS 


14.43 


20 





-350 


16 


13.46 



(0.09) 
0.10 
0.08) 



13.31 (0.09) 
13.74 0.10 
13.76 (0.09) 
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Table III. Table of optimal g 2 /4:ir and x 2 values obtained from x 2 maps. Solution A 
has only the one-pion exchange (OPE) contribution for J>6. Solution B has only the OPE 
contribution for J>5. Solution C is the same as A, with a plane-wave Born approximation 



for the OPE. The region below 400 MeV contains 2170 pp and 3532 np data. The region 
below 1600 MeV contains 12838 pp and 11171 np data. 



Solution 


g 2 /4n 
(pp data only) 


x 2 


9 2 /4n 
(np data only) 


x 2 


g 2 /4n 
(pp and np data) 


x 2 


A(0-400 MeV) 


13.61(0.09) 


3034 


14.16(0.12) 


4533 


13.80(0.07) 


7582 


B(0-400 MeV) 


13.70(0.08) 


3055 


14.07(0.11) 


4555 


13.83(0.07) 


7617 


C(0-400 MeV) 


14.42(0.08) 


3046 


14.51(0.11) 


4554 


14.38(0.06) 


7602 


(0-1600 MeV) 


13.67(0.06) 


22030 


13.42(0.04) 


17625 


13.51(0.03) 


39668 



12 



